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Abstract 

We examine the apphcabihty of perturbative QCD to B meson decays 
into D mesons. We find that the perturbative QCD formahsm, which in- 
cludes Sudakov effects at intermediate energy scales, is applicable to the 
semi-leptonic decay B Dlu, when the D meson recoils fast. Following 
this conclusion, we analyze the two-body non-leptonic decays B Dtt and 
B DDg. By comparing our predictions with experimental data, we extract 
the matrix element \Vch\ = 0.044. 



1. Introduction 



Recently, perturbative QCD (PQCD) has been proposed to be an alter- 
native theory for the study of B meson decays 0, |], which complements 
the approach based on the heavy quark effective theory (HQET) Q and the 
Bauer-Stech-Wirbel method 0. The point is to include Sudakov effects [|], 
which arise from the all-order summation of large radiative corrections in the 
processes. It has been shown that these effects, suppressing contributions due 
to soft gluon exchange, improve the applicability of PQCD to exclusive pro- 
cesses around the energy scale of few GeV . The heavy b quark possesses a 
large mass scale located in the range of applicability . The Sudakov factor 
for the heavy-to-light transition B —>■ ttIu has been derived in [Q, and the 
perturbative evaluation of the associated differential decay rate is found to 
be reliable for the pion energy fraction above 0.3. 

In this paper we shall investigate the applicability of the above PQCD 
formalism to heavy-to-heavy transitions, concentrating on the semi-leptonic 
decay B Dlu. Heavy quark symmetry ^ has been employed in the 
analysis of this deacy [0, whose amplitude is written as 



A{P^,P2) = ^V,,u^,{l-^,)l{D{P,)\c^%\B{P,)) , (1) 

where Gp = 10~^ GeV~^ is the Fermi coupling constant, and Pi (P2) is the 
B {D) meson momentum. The transition matrix element = (D|c7^6|i?) 
can be expressed in terms of a universal form factor ^ in the heavy meson 

limit i, 

iVP^ = y/memD^ivi ■ V2) {vi + V2T (2) 

with ruB {friu) the B (D) meson mass. The velocities Vi and V2 are defined 
by the relations Pi = msVi and P2 = m£,V2-, respectively. The form factor 
^, called the Isgur-Wise (IW) function, depends only on the velocity transfer 
Vi ■ V2, and is normalized to unity at zero recoil Vi ■ f 2 = 1 in the limits m^, 
mn ^ 00 0. 

The IW function has been usually regarded as sensitive to long-distance 
effects, and can not be calculated in perturbation theory. For the behavior 
of ^ above zero recoil, there is only the model estimation from the overlap 
integrals of heavy meson wave functions In this work we shall show 



that PQCD can give reliable predictions to ^ in the region with large Vi ■ V2, 
where the heavy quark symmetry can not provide any information of ^. We 
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argue that the IW function is dominated by soft contribution in the slow D 
meson hmit, at which the heavy meson wave functions strongly overlap, and 
factorization theorems do not hold. However, when the D meson recoils fast, 
carrying the energy much greater than mu, the case is then similar to the 
5 — i> TT decays, and PQCD is expected to be applicable |jT2| . 

The above conclusion then indicates that two-body non-leptonic decays 
such as B ^ Dn and B —>■ DDg can be analyzed reliably in the PQCD 
formalism. The B ^ D decays have been studied |jl2|, |T^ based on the 
exclusive PQCD theory developed by Lepage and Brodsky However, 
these analyses are lack of quantitative justification for the perturbative cal- 
culation, and are highly sensitive to the variation of the heavy meson wave 
functions. Our predictions for the branching ratios of these decay processess 
are comparable with those from the standard PQCD in [13[, and lead to 
the value 0.044 for the Cabibbo-Kobayashi-Maskawa matrix element \Vcb\ by 
combining with experimental data |T^. On one hand, we derive the behavior 
of the IW function near the high end of vi ■ f 2. On the other hand, the con- 
sistency of the extracted \Vcb\ with its currently accepted value justifies the 
application of our PQCD formalism to the semi-leptonic decays B —>■ nlu 
and i? — > TTTT 0]. 

A model-independent extraction of the matrix element \Vcb\ has been 
obtained from the semi-leptonic decay B D*lv in the framework of HQET 
[ p!0|l . The value of \Vcb\ was read off by extrapolating the experimental data 



to the zero-recoil point, at which the IW function is known to be equal to 
unity. In the present work, however, we must extract \Vcb\ by studying the 
behavior of the IW function at the opposite end of the velocity transfer, for 
which the PQCD analysis is reliable. Hence, the two-body decays are good 
candidates. Another possible method of extracting \Vcb\ has been proposed 
in |T^ , in which a sum rule for the relevant structure function of the inclusive 
semi-leptonic decay b ^ c was considered. 

In section 2 we derive the factorization formulas for the form factors in- 
volved in the B —>■ Dlv decay, including the resummation of large radiative 
corrections to this transition. Numerical analysis is shown in section 3, along 
with the behavior of the IW function at large velocity transfer. The com- 
parision of our predictions for the decays B Dn and B DDg with 
experimental data is also made. Section 4 contains the conclusions. 
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2. Factorization 



In this section we develop the factorization formula for the B Dlu 
decay. The lowest-order factorization for the transition matrix element 
is shown in fig. 1, in which the bubbles represent the B and D mesons, 
and the symbol x represents the electroweak vertex where the lepton pair 
emerges. The b quark, denoted by a bold line, and its accompanying light 
quark carry the momenta Pi — ki and ki, respectively, which satisfy the on- 
shell conditions (Pi — ki)^ ^ ml and kf ^ 0, rrih being the b quark mass. 
We shall work in the rest frame of the B meson such that the nonvanishing 
components of Pi are Pi^ = Pi~ = ttib/V^- contains small amount of 
transverse components ki^-, and its minus component defines the momentum 
fraction Xi = k^ / P{ . The assignment of the momenta for the D meson 
is similar, but with /ci, rrii, and xi replaced by /c2, 'mc and X2 = k2 / P2 , 
respectively, rric being the c quark mass. 

The expressions for the components of P2 are more complicated. At zero 
recoil the D meson sits at rest with the B meson, and we have P2 oc Pi. 
When the D meson takes the maximum energy, it moves fast and P2^ is 
much greater than P2"". To show the relation between P2^ and P2~, it is 
most convenient to express them in terms of the velocity transfer rj = vi ■ V2- 
Solving the equations Pi ■ P2 = rjmBTnj:, and P2 = mj^, we obtain 



7] + - 1 

71" 



P/ = ruD 



P2 = mo . (3) 

The upper bound of r/, corresponding to the maximum recoil of the D meson, 
is equal to r/max = {rnB/rriD + ijid IruB) l'2.. It is easy to check from eq. (|^) 
that P2'*' = P2~ = m£,/\/2 as r] takes the minimum value 1, and P2 / P2 = 
m\/m\, > 1 at ?7 = ?7max- 

We then consider higher-order corrections to the basic factorization pic- 
ture. As analyzed before |]^, ^ |^, these corrections produce large collinear 
logarithms, when the loop momentum is parallel to that of a light quark, or 
large soft logarithms, when the loop momentum is much smaller than the 
mass scale involved in the processes. The two types of large corrections may 
combine to give double logarithms. It has been shown that the double loga- 
rithms come from two-particle reducible diagrams in physical (axial) gauge. 
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whose contributions are dominated by collinear enhancements for fast hght 
mesons, and are dominated by soft enhancements for heavy mesons at rest 
[0]. Therefore, they can be absorbed into the corresponding wave functions, 
which involve similar dynamics. The all-order summation of the double log- 
arithms in light meson wave functions, such as a pion, has been performed in 
[0]. The resummation technique Jl^ has been extended to the case of heavy 
mesons in [|]. Combining the above results, we have derived the Sudakov 
factor for the heavy-to-light transition B ttIu [|]. 

The analysis of Sudakov corrections to fig. 1 is more complicated com- 
pared to that of the decay B tcIu. Due to the dominance of soft contri- 
bution near the low end of rj, we concentrate only on the large t] region. In 
this region radiative corrections on the D meson side involve three scales, 
P2' rrif) ^ k2T- Note that all the previous studies of resummation in- 
volve only two scales, for example, and fc-r in the pion case, and 
and kT in the B meson case. The three scales produce various large loga- 
rithms of P2 /k2Ti P2 I^D and mj:,/k2Ti which complicate their organiza- 
tion. As a noive approximation, we keep only the largest one proportional 
to \n.{P2 /k2T)- The neglect of those logarithms containing is equivalent 
to the neglect of P2 ^ P2 in the analysis of radiative corrections to the D 
meson wave function. The D meson is then regarded as light in the large r] 
region, and the corresponding Sudakov factor for the decay B Dlu can be 
approximated by that for the heavy-to-light transitions. 

The factorization formula for in the transverse configuration space, 
with radiative corrections taken into account, is written as 

X H''{xi, X2, bi, b2, niB, tjid, fi) Vb{xi, bi. Pi, /i) , (4) 

in which both the B and D meson wave functions, Vb and Vd-, contain 
the evolution from the resummation of double logarithms performed in axial 
gauge. We have introduced the conjugate variable hi (62) to denote the 
separation between the two valence quarks of the B (D) meson. We shall 
employ the approximation mb ^ tub = 5.28 GeV and rric ~ rriD = 1.87 GeV 
in eq. H'^ is the Fourier transform of the hard scattering amplitude 
to b space, fi is the renormalization and factorization scale. 

Note that in the evaluation of H'^ we neglect those terms proportional to 
ki and A;^ in the hard scattering amplitude following the kinematic ordering 
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~ ^2 ^1 ~ ^2^, which is vahd in the large rj region. For example, the 
gluon propagator in the lowest-order diagram is written as 

(5) 



where serves as the infrared cutoff of the Sudakov corrections. Once the 
approximation is made, the and /c^ dependences, appearing only in the 
B and D meson wave functions, respectively, are integrated to give eq. (Q). 

As stated above, near the high end of i] the Sudakov factor exp(— 5") for 
the decay B Dlu, which groups the large logarithmic corrections in Vb, 
Vd and H^, can be approximated by that for the heavy-to-light transition 
derived in with the exponent S given by 



2 

X=X2, 1—X2 



S{xi,bi,mB,mD) = s{xi,bi,P^)+ ^ s{x,b2,P2 

(6) 



1 



1^ ln(t/A) ^ ln(t/A) 



ln(6iA) -ln(62A) 



where t is the largest mass scale associated with the hard gluon, and will 
be specified later. The first term in eq. comes from the resummation of 
redicible corrections to the heavy meson wave function The value of A = 
Aqcd will be set to 100 MeV below. The complete expression for the factor 
s{x,b,Q), including the leading and next-to-leading logarithms, is exhibited 
in Appendix. It is observed that exp(— S") decreases quickly in the large bi 
region and vanishes as bi > 1/A. Therefore, the long-distance contributions 
are suppressed, and the perturbative calculation becomes relatively reliable. 

One may wonder whether the resummation of large radiative corrections 
can improve the applicability of PQCD near the low end of r^. If we recog- 
nize that the D meson is regarded as a heavy meson in this region, and is 
dominated by similar dynamics to that of the B meson, the Sudakov factor 
for the decay B — >■ Dlu can be taken as the combination of the expressions 
for heavy mesons @] at two different mass scales, rriB and tjid- The Sudakov 
exponent S is then written as 



S{xi,bi,mB,mD) = s(xi, 6i, ) s(x2, 62, 

1 



ln(t/A) ^ ln(t/A) 



ln(6iA) -ln(62A) 



(7) 
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Obviously, it is not expected that our perturbative analysis with the above 
Sudakov suppression becomes self-consistent. The virtuality of the hard 
gluon in fig. 1 diminishes as Xi and X2 are both small, which leads to a 
large running coupling constant as- However, this non-perturbative region 
is not strongly suppressed by the Sudakov factor in eq. (^. It is the extra 
exponent s(l — X2, b2, P^) in eq. (H) that can provide necessary suppression 
in the small X2, or large 1 — X2, region. 

Having factorized all the large logarithms into the Sudakov factor, we can 
then compute the hard scattering amplitude H^^ of the B Dlv decay to 
lowest order of a^. From fig. la we have 



tr 



X 



75(/p2 + ^d) f, /Pi- M + rriB a {A + "^i?)75 
Vm ^ {Pi-k2f-ml^ 

ih - k2Y 

IGnasCplmBrnD — X2(imy\ 
[xiX2CmBmD + (kiT - k2r)^] [x2C"^Bmz) + k^^' 
^ lQT[asCF[m% + ^2(2^3^,0] 

[xiX2(mBmr) + (ki^ - k2Tf][x2(mBmo + kj^ 



■-'2 ) 



with 



Ci 

C2 



+ 



V 



2 ■ 2^/^i^^ 

2r/2 -27]-l 



(9) 



The factors (/Pi + raB)^hl \P^^c and 75(/P2 + raD^j ^plW^ come from the 
matrix structures of the B and D meson wave functions, respectively. Cp = 
4/3 is the color factor, and Nc the number of colors. Similarly, from fig. lb 
we get 



tr 



X 



ih - k2f 
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ITl 



[xiX2CmBmD + (kiT - kar)^] [a;iC?TiB"^D + k' 
IQti a sCp[mBmj:) — XiC^imW 
[xiX2CmBmD + (kiT - \<.2TY][xiCmBmD + k' 



-A 



ITJ 



-Pa'^ . (10) 



Note that H^^'^ can be obtained from by exchanging the variables as- 
sociated with the B and D mesons. This permutation symmetry has been 
displayed manifestly in fig. 1. 

Performing the Fourier transform of eqs. (|^) and (|TDp to get and 
substituting them into (^), we obtain the factorization formula for = 
fiPi + /2-P2'! where the form factors /i and /a are given by 

/i = IQttCf / dxidx2 / bidbib2db2 (f)B{xi,bi)(j)D{x2,b2) 
Jo Jo 

X [{mBrriD - X2Cimjj)h{xi, X2, &i, ^2) 

+ (m|, + Xi(2mBmD)h{x2, Xi, 62, &i)] 
X exp[-S{xi, bi, rriB, mo)] , (11) 



and 



/2 = IQttCf / dxidx2 / bidbib2db2(l)Bixi,bi)(j)nix2,b2) 
Jo Jo 

X [(m| + X2C2mBmD)h{xi,X2, &i, ^2) 

+ {mBmD - XiCim|)/i(x2, Xi, 62, &i)] 
X exp[-S'(a;i, 6i, m^, mo)] , (12) 



respectively, with 



h{xi,X2,bi,b2) = as{t)Ko [JxiX2CmBmDb 



X 



9{bi - b2)Ko J X2CmBmDbi Iq J X2CmBmDh 



+6{b2 - bi)Ko \ \l X2C,mBmDb2j lo \ \l X2CmBmDbi 

(13) 

The wave function (pB {4'd) comes from Vb {Vd) in eq. @) with the evolution 
in Pi {P2~), which is the result of the resummation of reducible corrections, 
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grouped into the Sudakov factor. The argument 6 in 0b and (pD denotes the 
intrinsic transverse momentum dependence of the wave functions [1^, which 



is a non-perturbative object, and can not be handled in perturbation theory. 
Kq and Jo are the modified Bessel functions of order zero. We choose t as 
the largest scale associated with the hard gluon, 



t = max K/xia;2C'7J.B'7iD, 1/62 • (14) 



3. Numerical Results 



Before evaluating /j, we compare our formulas with those derived in the 
framework of standard PQCD |]12|, |l^, where the k-p dependence in the hard 
scattering amphtude is neglected, and the heavy meson wave functions, with 
ky integrated, take the simple form of the (5-function (the so-called peaking 
approximation) , 



B 



2VS 



6{x-xb), 4>d{x) 



D 



2^3 



6{x — xd) 



(15) 



Here Jb = 0.12 GeV and = 0.14 GeV are the decay constants of the B 
and D mesons [0, respectively. Eqs. (|1T]) and (p!^ ) are then reduced to the 
standard factorization formulas without b integrations, which lead to 



fi 



f2 



4 

-T^CpasfsfD 
4 

-T^CpasfsfD 



rriBmr, - XoCim]^ _^ "^z? + XBC2mBmD 



XBX^C^m^m'j^ 
m% + XDC2mBmD 
XB^jjC^m'^mj^ 



+ 



x^XoC'^rn'^m'jj 
XBCim% 



niBniD 



XB^oC^fTlB'^D 



(16) 



It is apparent that the above expressions are very sensitive to the values 
oi xb and x^, and the coupling constant as must be regarded as a free pa- 
rameter. We consider the non-leptonic decay B Dn, which corresponds 
to the case of maximum recoil here with rj = rimax = 1-59. Setting as = 0.4, 

= 1.3, which gives a 
However, if slightly 



xb = 0.07 and xd = 0.2 as in [|T3l, we obtain /i + /2 



branching ratio comparable with experimental data |T 
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different values like xb = 0.07 and xd = 0.15 were inserted, the branching 
ratio becomes 3 times larger. On the other hand, simply setting as to a 
constant makes the justification of the perturbative calculation unavailable. 
Compared to the standard PQCD approach, our modified perturbative ex- 
pressions are less sensitive to the profile change of the wave functions due to 
the inclusion of in the hard scattering amplitude, which moderates the 
divergences from small xb and xjj. Substituting eq. (^) into (PD and ([T^) , 
and performimg the integrations over bi and 62, we find that the latter set 
of Xb and Xd leads to a branching ratio only 50% larger than that from the 
former set. 

We adopt the following model for the B meson wave function [^ , 

-2 



$ij(x,k2 



B 



Cb + 



m 



B 



+ 



x] 



1 — X x{l 

The constants Nb and Cb are determined by the normalizations 

fB 



(17) 



which give Nb 
defined by 

(pB{x,b) ■■ 



/^'dx/d^k.M^,k.)-^^, 

|^'dx|d2kT[$B(x,kT)F = l, 

0.923 GeV=^ and Cb = -27.877255 GeV^. 



d^kT^B{x,kT)e'^^-^ 

nNBbx'^il - x f 



(18) 
')B is then 



m%x + Cbx{1 — x) 



Ki Umlx + CBx{l-x)b] . (19) 



It is observed that (pB peaks at s ^ 0, and decreases monotonically with x 
for a fixed b, signifying the soft dynamics involved in the rest B meson. 

If assuming the similar model for the D meson wave function with in 
eq. ([T7| ) replaced by m/j straightforwardly. 



<pD{x,b) 



nNobx^n 



X) 



mj-fX + Cr)x{l — x) 



Ki \Jm\)X + Cdx{1 — x)b 



(20) 



we obtain the constants Nd = 0.136 GeV^ and Cd = -3.495345 GeV^. 
The resulting wave function (po also peaks at small x ~ 0.01 for a fixed 
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b. However, the QCD sum rule analysis in |^ has shown that the average 



momentum fraction of the light valence quark in a fast D meson is roughly 
0.2. To be consistent with this observation, we employ eq. ( PD| ) but with Cd 
determined by the requirement that (p^ takes the maximum value at x ~ 0.2 
for b ^ 0. We then have Cd = -2.9 GeV^, along with A^^ = 0.240 GeV^ 
from the normalization / da;0£)(x, 0) = /d/(2v^). 

Results of /i and /2 derived from eqs. ([TID and (|12D, respectively, with 
bi and 62 integrated up to the same cutoff be are shown in fig. 2. We find 
that at ?7 = 1.30 approximately 55% of the contribution to fi comes from 
the region with as{l/bc) < 1, or equivalently, be < 0.5/A. The percentage of 
perturbative contribution increases with rj, and for rj above 1.39, more than 
60% of the full contribution is accumulated in this region. It implies that our 
PQCD analysis of the decay B Dlv in the range of r/ > 1.39 is relatively 
reliable, since perturbative contribution dominates [0]. It is also found that 
the self-consistency of the perturbation theory becomes worse quickly for 
T] < 1.3 as expected. 

Based on the above conclusion, we are led to consider the two-body non- 
leptonic decays such as -B ^ Dtt and B DDg, which can be best described 
by our PQCD formalism. The decay rate of the specific mode B^ — > D^7i~ 
is given by 



-GlK,\WfM (1-^1 1/1 + /2P , (21) 



647r ' rnly 

which is derived from the amplitude 

A=^V^,Veb{7r\j,{l-j,)\0){D\crb\B) (22) 

with the PCAC relation (7r(P)|7^(l - 75)10) = iV^f^Pf, inserted, = 93 
MeV being the pion decay constant. Eq. (p2D is achieved following the con- 



clusion in ||2^ that the non-factorizable ly-exchange contribution is negligi- 
ble. The value of /i -|- /2 in this case can be easily read off from the curves 
corresponding to r] = 1.59 in fig. 2, which is equal to 1.44. Substituting 
the matrix element = 0.974, we obtain F = 8.4 x 10"-^^|Vcbp GeV, or 
equivalently, the branching ratio B{B^ D^n^) = 1.65\Vcb\'^ from the total 
width (0.51 ±0.02) x 10~^ MeV of the B^ meson ||23l. Comparing with exper- 
imental data B{B° D+ TT ) = 3.2 X 10 ^, we extract the matrix element 
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\Vcb\ = 0.044, consistent with currently accepted value ||23[. Similarly, the 
decay rate for the mode 5° — > D^D~ is given by 



327r 



m 



Vmax ^ 



B 



X 



'2 2, 2\Y-|/2 2 2 \ _[ 



(23) 



with the matrix element \Vcs\ = 1-0, the decay constant of the Dg meson 
ff)^ = 0.16 GeV [jl9l, and the meson mass niD^ = 1.97 GeV. In this case we 



have the maximum velocity transfer r]'^^^ = {m% + m'jj — mj^ J / (217131710) 



1.39, for which the corresponding values /i 



OB -r- "I'D 

0.47 and /s 



1.32 are read 



off from fig. 2. Eq. (^) then gives the decay rate F = 2.7 x 10 ^^|Vcfep, or 
the branching ratio B{B^ D^D^) = 5.3\Vcb\'^. Experimental data show 



B{B^ D+D-) = 9.9 X 10~^ from which we extract = 0.043, close to 
that obtained from the decay B^ D~^tt^. 

Due to the consistency of our predictions with experimental data, we can 
explore the behavior of the IW function near the high end of rj reliably. For 
finite and m^, eq. (|]) is modified to 

> in the heavy meson limit. A simple manipulation 



where C,+ ~^ ^ and ^ 
gives the relations 




(25) 



The dependence of ^+ and on rj is shown in fig. 3, which exhibits a falloff 
and an increase with i], respectively. The magnitude of ^_ indeed diminishes 
as stated above. A model calculation of the IW function in terms of the 



overlap integrals of the heavy meson wave functions has been performed [11 
which leads to 



exp 



-(2p2-l 



(26) 



1] + 1 [ 1] + l_ 

with the parameter p ^ 1. The behavior of C,Tn is also shown in fig. 3. It is 
observed that our predictions for ^+ are close to at large r], and begin to 
deviate from as 77 < 1.39. The match confirms the applicability of PQCD 
to heavy meson decays in the large recoil region. 
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At last, the differential decay rate for the specific mode D^l v 

with vanishing lepton masses is given by 

^ ^ \V.,fR{^) = \V.,\'^,GlmBmW - ifVi + • (27) 

Substituting the results of fi into eq. (|27|), we derive the behavior of R{ri) 
for rj > 1.39 as in fig. 4, which shows an increase with rj. Once experimental 
data for the spectrum of the decay —>■ D~^l~u are available, the matrix 
element \Vcb\ can also be extracted from the curve in fig. 4. 



4. Conclusions 

In this paper we have applied the PQCD formalism to the semi-leptonic 
decay B Dlu, and found that the perturbative calculation is reliable for 
the velocity transfer above 1.4. The point is to include the resummation of 
large radiative corrections in the process, which improves the applicability of 
PQCD. The intrinsic transverse momentum dependence also plays an essen- 
tial role in the calculation. We emphasize that our analysis does not involve 
any phenomenological parameter, and is insensitive to the profile change of 
the wave functions. The perturbative calculation is justified by considering 
the magnitude of the running coupling constant, which defines the region 
where the perturbation theory is reliable. 

Our predictions are satisfying in the sense that they match the model 
estimation of the IW function at the high end of velocity transfer, and the 
values 0.044 and 0.043 for the matrix element \Vcb\ are extracted from the 
decays B Dn and B DDs, respectively. The results presented in this 
work confirm our perturbative analysis of the decay B tiIu in 0, which is 
important for the extraction of the matrix element \Vub\- 

I thank F.T. Cheng, T. Huang, G. Sterman and H.L. Yu for helpful 
discussions. This work was supported by the National Science Council of 
R.O.C. under Grant No. NSC84-2112-M194-006. 
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Appendix 



In this appendix we show the derivation of the exponent s{x, b, Q) in 
eq. §). We start with eq. (5.42) in ref. 



Q d/i 

l/b fi 



In 1-^1 A{g{^i)) + B{g{^i)) 



s{x, b, Q) 

in which the factors A{g) and B{g) are expanded as 



(28) 



Big) 



hi 

3 71 



(29) 



in order to take into account the next-to-leading logarithms. The running 
coupling constant a, is written as 



a,(/i) 



1 



/32lnln(^VA2^ 



TT /?iln(/iVA2) /3f ln2(/i2/A2) 

The above coefficients /3i and A^^^ are 

33 - 2n^ 153 -19n; 

Pi — — tt; — 5 P2 — - 



(30) 



12 



^(1) 



^(2) 



24 

67 10 



9 



27 



+ In 



(31) 



where tt,/ = 4 is the number of quark flavors, and 7 is the Euler constant. 
Performing the integration in eq. (P^, we obtain s, which is given in terms 
of the variables 



g = ln(a;Q/A), 6 = ln(l/6A) 



(32) 



by 



AW 



AW 



AW 



aw aw 

4/3? 4/3i 



In 



,27-1 ■ 



In 



4/3? \ 'b 

( q 



4 - 1 



13 



ln(2g) 



In (26) + 1 



q b 

\n\2q) ~ln\2b) 



In 



,27-r 



ln(2g) + 1 ln(26) + 1 



21n(2g) + 3 21n(26)+3 



2q 



21n(26) 



18 ln^(2g) + 30 ln(2g) + 19 18 ln2(26) + 30 ln(26) + 19 



A^^^f3lq-'b 
432/3f S3 



62 



91n^(26) + 61n(26) +2 



(33) 



The previous studies involving the Sudakov logarithms pick up only the first 
six terms in eq. (|33|) , which are more important than the remaining ones in 
large Q region. Note that the coefficients of the fifth and sixth terms are 
different from those in refs. @, 0. It can be easily checked that with these 
corrections the results for the pion form factor in |^ are reducend only by 
few percents. An explicit examination on the form factors fi in B ^ D 
decays shows that the partial expression including only the first six terms 
give predictions smaller than those from the full expression by less than 
5%. Hence, for simplicity this partial expression is substituted into ([TID and 
(p!2D. Note that s is defined for q > b, and is set to zero for q < b. As a 
similar treatment, the complete Sudakov factor exp(— S*) is set to unity, if 
exp(— S) > 1, during the numerical analysis. 
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Figure Captions 

Fig. 1. Lowest-order diagrams of tlie decay B —>■ Dlv. 

Fig. 2. Dependence of (a) /i and of (b) /2 on the cutoff he for (1) 77 = 1.3, 
(2) 7] = 1.39, and (3) r] = 1.59. 

Fig. 3. Dependence of (a) ^+ and of (b) on rj derived from our PQCD 
formalism. The dependence of C,m on rj from the model calculation in []Tll 
(dashed line) is also shown. 

Fig. 4. Dependence of R{r]) on rj derived from the PQCD analysis. 
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